Non-thermal equilibration of a one-dimensional Fermi gas 
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Equilibration of an isolated Fermi gas in one spatial dimension after an interaction quench is 
studied. Evaluating Kadanoff-Beym dynamic equations for correlation functions obtained from the 
two-particle-irreducible effective action in nonperturbative approximation, the gas is seen to evolve 
to states characterized by thermal as well as nonthermal momentum distributions, depending on the 
assumed initial conditions. For total energies near the Fermi temperature, stationary power laws 
emerge for the high-momentum tails while at lower momenta the distributions are of Fermi-Dirac 
type. The relation found between fluctuations and dissipation exhibits nonthermal final states. 
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The dynamics of interacting quantum many-body sys- 
tems far from equilibrium remains a challenge both for 
theory and experiment. One of the seemingly simple 
questions concerns the long-time evolution of an isolated 
such system: Given that it starts in some deliberately 
chosen initial non-equilibrium configuration, will the in- 
teractions between particles force the system to equili- 
brate, and if yes, to what kind of state? While for clas- 
sical integrable systems such problems have been stud- 
ied at length since the pioneering work of Fermi, Pasta, 
and Ulam [l[ , the evolution of quantum systems remains 
much more elusive. Does a closed and finite but suffi- 
ciently large system with a ground state show, at least 
approximately, equilibration to a thermal or to some al- 
ternative quasi-stationary state before its unitary evolu- 
tion causes revivals 0, and if so, on what time scale does 
this happen? In recent years, new experimental tech- 
niques have been developed for ultracold atomic gases 
and solid-state physics that allow to precisely study quan- 
tum many-body dynamics. Experiments have focused on 
the question of equilibration of a one-dimensional (ID) 
quantum gas with quadratic dispersion in which isolated 
binary collisions of the particles would not change their 
momenta owing to the restrictions of simultaneous mo- 
mentum and energy conservation [3, |3| • 

In this article, we study the long-time evolution 
and equilibration of a ID Fermi gas containing two 
spin components a € {tii} that mutually interact 
through local repulsive s-wave collisions described by 
the Hamiltonian H = J dx['ifl^{x){~d'^/2m)'i>a{x) + 

igo.f)/2)^i{x)^l{x)'fi3{x)'i/o.{x)], with go^p = (1 - (5a/3) 
AiTaiD/m, fliD being the ID scattering length, x ~ 
(xo,x) includes time and space, and the complex 
fermionic fields obey canonical anticommutation rela- 
tions {^'a(x, t), \E'^(y, <)} = (5q^(5(x — y). This model is 
integrable and has as many conserved quantities as there 
are degrees of freedom [5| such that it is expected not 
to thermalize in general, i.e., not to approach a grand- 



canonical ensemble [6|-|8[. Its low-energy properties can 
be approximated by a Tomonaga-Luttinger liquid (TLL) 
showing, e.g., a non-Fcrmi-Dirac power-law momen- 
tum distribution around the Fermi surface. The long- 
time evolution of the TLL model, following an interaction 
quench, was found in [l^l to approach a generalized Gibbs 
ensemble @ determined, besides energy, by the values 
of conserved quasiparticle numbers. Such non-thermal 
stationary states were also found in other models, both 
integrable and non- integrable ones 11 -[l4 1 . 

The results of evaluating quantum-field theoretical 
methods wc present here imply that a ID Fermi gas 
defined by the full non-quadratic Hamiltonian H thcr- 
malizes for sufficiently large initial energies while it can 
approach a non-thermal state at lower energies, with ul- 
traviolet power-law tails in the momentum distribution 
away from the Fermi edge. We consider the dynami- 
cal evolution of two-point correlation functions, including 
single-particle densities, as described by the Schwinger- 
Dyson or KadanofF-Beym equations 
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These equations have been used to describe thermal- 
ization of both relativistic and nonrelativistic quan- 
tum gases 15, [l^. The statistical (F) and spec- 



tral (p) components of the connected two-point Green 
function are defined as Gab{x,y) = {T^l_{x)'i'b{y)) = 
Fab{x,y) - (i/2)pab{x,y) sgn(a;o - yo) where T denotes 
time-ordering, sgn(i) is the sign function evaluating to 
+1(-1) for t > {t < 0). The field index a = {a,i) 
combines the spin index a and the index i numbering 
the complex fields ^a,i{x) = '^aix) and ^'0,2(2;) = 
vl'jj,(x) into a single index. The mean-field energy is 
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FIG. 1: (Color online) Diagrams contributing in leading order 
(LO) and next-to-leading order (NLO) of the 1/M expansion 
to /2 in the two-particle-irreducible effective action. Bare 
vertices are depicted as wiggly lines at the ends of which it is 
summed over spin and field indices a. Vertices evaluate the 
incoming (blue) lines denoting G at a single space-time point 
over which it is integrated. Statistical factors are omitted. 

Mab{x) = Sab[-d^/ {2m) - gajFccix, x) /2]+ gapFabix, x), 
where sums over spin 7 =tj4- ic = 1,2, with 

c (7, ic) are implied. The right-hand-sides of Eqs. ^ 
contain the non-Markovian memory integrals over the 
correlation functions with J^^ ^ = J^^ dzQ J dz. The sta- 
tistical and spectral components of the "self-energy" 
Sab{x,y) = S^i,{x,y) - {i/2)S^^^{x,y) sgn(x-o - yo) con- 
tain all the information about the correlations in the 
system due to interactions between particles during the 
past evolution. Using a variational principle, we de- 
rive the above dynamic equations from the two-particle- 
irrcducible (2PI) effective action r[G], a functional of 
G, see, e.g., Refs. [H, [l^. One has Sab{x,y) = 
—2iSr2[G]/SGba{y,x) of the beyond-one-loop part /2 of 
this effective action, and Sab{x,y) contains all diagrams 
beyond mean- field order, S = T^^) + S, while the local 
part IJ^°^ {x) contributes to M. 

Different truncation schemes are available, which we 
found for our cases to yield qualitatively and quantita- 
tively close results. Fig. [1] shows the leading-order (LO) 
and next-to-leading-ordcr (NLO) contributions of the ex- 
pansion in 1/Af we take into account, where TV is the 
number of different spin states [l3| • In order to check the 
dependence on the truncation, we compared this with a 
truncation including all order-t/^ loop diagrams. 

We solve Eqs. ([T]), together with the integral equations 
determining S ' ,m one spatial dimension in a finite-size 
box with periodic boundary conditions. The initial state 
is assumed to be entirely defined by the Green function G. 
Working in Fourier space, the initial state is encoded into 
the Go&(xo,yo;p) = / (ircxp[-ipr]Ga6(xo, r; 2/0, 0) at the 
initial time xq = yo = Iq. Specifically, we fix the single- 
particle momentum distribution (cr^ = Pauli 3-matrix) 

n„(t,p) = [l + crfjF(^c.,j){a.t)it,t;p)]/2, (2) 

obtained through the definition of the statistical cor- 
relator Fab{x,y) = ([*l(a;),*b(2/)])/2, at time t = to 
while the spectral function pab{x,y) = i({5'Jj(a;), 5* (,(2/)}) 
is fixed for xq ~ yo by the fermionic anticommutator, 
Pabit-,t;p) = i6ab- Furthermore, we choose the coher- 



ence between different spins as well as the pair corre- 
lation function to vanish, (\D'J,(<o,x)^'/3(to,x)) — 0, for 
a ^ 13, and (*Q(to, x)*/3(to, x)) = 0. This implies 
F{a,i){p,j){ttt\P) = 0, for * 7^ j- We restrict ourselves, 
here, to the special case of spin-balanced systems. 

A prominent strength of the 2PI approach to far-from- 
equilibrium many-body quantum dynamics is that the 
equations of motion are derived from an action functional 
and therefore, for a closed system, automatically conserve 
total energy and further crucial quantities like the local 
particle current-density vector [16|. As a consequence, 
the dynamic equations considered here allow to study the 
long-time dynamics and equilibration of a system starting 
far away from equilibrium. 

We first present results where the chosen initial con- 
ditions allow the single-particle momentum distribution, 
within the range of momenta considered, to thermally 
equilibrate. We find that the values of the conserved 
quantities present in the initial state determine the fi- 
nal state while all other information about the details of 
the initial state are lost during the evolution. Fig. [2^ 
shows two different initial distributions with the same 
total particle number and energy from each of which the 
system approaches the same final state. The coupling 
strength is chosen as 7 = mg/niYi = 4. Off-energy-shcll 
two-body scattering in the many-body background al- 
lows the momenta to be redistributed despite asymptotic 
conservation of momenta in isolated two-to-two scatter- 
ing. Fig. [2)3 shows the occupation numbers as a func- 
tion of time t for six different momentum modes. The 
evolution is characterized by a short-time scale of de- 
phasing depending on the width of the momenta present 
initially and a long-time scale determined by scattering. 
A comparison with classical hard-sphere scattering gives 
that the equilibration takes on the order of 3 to 4 colli- 
sion events. We note that, for typical line densities of 
riiD — 10^ cm^^, the thermalization time of a Li gas 
results as Tth — 0.25 ^s. Taking into account that the 
thermalization rate scales approximately as (oid^^id)^ 
this implies thermalization times as large as seconds for 
a weakly interacting gas with the same niD. 

To exhibit the Fermi-Dirac character of the final mo- 
mentum distribution, we show, in Fig. [2j:, the evolution 
of the inverse-slope function /(p) = ln[l/ ^.-[-(p) ~ 1] i^i run 
B. For a Fermi-Dirac distribution n^(p) ~ ??fd(i^p — /i), 
nFD(w) = [exp(/3a;)-|-l]~^, this function is represented by 
a straight line with slope /3 = l/fc^T determined by the 
inverse temperature and ordinate-intercept — /?/-t related 
to the chemical potential p.. 

To investigate the character of the state approached 
at large times further, we compared the asymptotic mo- 
mentum distributions obtained for different initial ener- 
gies, for the same total particle number and interaction 
strength as before. Fig. [2^ shows the initial momen- 
tum distribution for run B and two further such distribu- 
tions C and D. Solving the equations of motion (IT|), we 
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(a) p/pf (b) nl^t/m (c) w(p)/a;_p 



FIG. 2: (color online) Equilibration of a one-dimensional Fermi gas starting from different initial momentum distributions, 
(a) Initial (t = 0) and final {t — 10 mn^^) distributions n^(t, |p|) — n_i{t, |p|) in runs A and B have the same total particle 
number and energy. In runs C and D, the particle number is the same while the energies are lower than in B. (b) Occupation 
numbers n-f{t, |p|) as a function of time t for momentum modes Pi = sin[j7r/A^s]/as, Ns — 128. (c) Inverse-slope function 
cr(a;p) = \n[l/nf{t,p) — 1] of the distribution nf{t,p) as a function of mode frequency iUp in units of the Fermi frequency ujf 
during the thermalization process, a reduces to a straight line for a Fermi-Dirac distribution n-t-(Lj) = {exp[/3(aj — /i)] + 




FIG. 3: (color online) (a) Inverse-slope functions a for runs B, C, and D at late times. The black dashed lines are Fermi-Dirac 
fits to the lowest 14 momentum modes from which the temperatures and chemical potentials are extracted, (b) Dependence 
of the chemical potential /i on the final temperature T — (fcs/?)~^, both extracted from a for various runs with the same total 
particle number but different energies, (c) Late-time momentum distributions, equivalent to those shown in panel (a). 



find that all runs reach a stationary momentum distribu- 
tion; however, for lower initial energies, they do no longer 
correspond to a Fermi-Dirac distribution over the entire 
range of momenta. As shown in Fig. [3^, the lower mo- 
mentum modes thermalize while the higher momenta in 
runs C and D remain overpopulated as compared to the 
exponential fall-off of the Fermi-Dirac distribution. The 
long-time approach to these distributions is qualitatively 
similar to the ones shown in Fig.l^b, i. e., the distributions 
appear frozen on a double-logarithmic scale. This leads 
us to conjecture that they do not undergo significant fur- 
ther changes over a long time. Note that the difference 
between the observed and thermal high-momentum occu- 
pation numbers for the respective temperatures is several 
orders of magnitude. 

Keeping the total particle number constant, a further 
reduction of the population at high momenta would not 
significantly alter the population of the lower momenta. 
Hence, we extract temperatures and chemical potentials 
from a fit of the lowest 14 momentum modes to a Fermi- 
Dirac distribution. The temperature dependences of the 
so found chemical potential /i is shown in Fig. [HJd. At 
high T, the chemical potential of the interacting gas (7 — 



4) converges to the value found for an ideal Fermi gas 
in thermal equilibrium on a finite-size discrete lattice, 
clearly distinct from the behavior in the thermodynamic 
limit. At low temperatures, we find significant deviations 
from the ideal- gas behavior as shown in the inset. 

In Fig. [3];, we depict the momentum dependence of 
the occupation numbers for runs C and D on a double- 
logarithmic scale. This shows that the overpopulation at 
high momenta is characterized by a power-law n{p) oc 
with K ~ 4.4. This exponent does not change when 
we include order-i?^ diagrams. We point out that, com- 
pared to the situation in a TLL, we find algebraic fall-off 
at high momenta away from the Fermi edge. 

This power-law tail in the momentum distribution 
suggests but does not automatically prove that the 
one-dimensional Fermi gas approaches a non-thermal 
state. To investigate this point further, we have a closer 
look at the fluctuation-dissipation relation between the 
Fourier transforms of the statistical correlation function, 
Fqq(Xo;i:j,p) = J dsexp{iu;s)FaaiXo + s/2, Xq - s/2; 
p), Fca{t,t';p) = {[^l{t,p),-^cit',p)])/2, and the spec- 
tral function pac.{t,t';p) = p), ^'^(t', p)}). Pro- 
vided that the many-body quantum state represents a 
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FIG. 4: (color online) Upper row: Spectral functions as a 
function of frequency at late time Xq = 18.9 n~[^m in runs C 
and D for five of the momentum modes pi, see legend in lower 
right panel. Lower row: Inverse-slope function of fractions / 
of the statistical F divided by the spectral function p at Xo — 
18.9 n^^m, for the same five momentum modes, see main text. 
Black lines indicate Fermi-Dirac distributions with /3 and fi 
as in Fig. [3^. In run D, the system does not thermalize. 



elude that the system thermalizes to a grand-canonical 
ensemble, with the eigenmodes of the strongly interacting 
system being superpositions of particles and holes. 

Run D, however, performed at even lower energy, 
shows that the system does in general not thermalize 
to a grand-canonical ensemble. As shown in the lower 
right panel of Fig. 21 the result of this run violates the 
fluctuation-dissipation theorem. Although the momen- 
tum overpopulation is again largely produced by the con- 
tributions from the Fermi sea, see Fig. [4] (upper right 
panel), also the fraction / shows a power-law tail ~ p^^. 
In conclusion, we find dephasing and equilibration of a 
ID Fermi gas with positive contact interactions, at suffi- 
ciently high total energies, to a thermal state satisfying 
the fluctuation-dissipation theorem for a grand-canonical 
ensemble, while at lower total energies non-thermal sta- 
tionary power-law distributions arc found. 

We thank C. Bodct, M. HoUand, S. Jochim, S. Kehrein, 
B. Nowak, J. M. Pawlowski, A. M. Rcy, and D. Sexty 
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Forschungsgemeinschaft, the Alliance Program of the 
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the German Academic Exchange Service (DAAD) is 
acknowledged. 



grand-canonical ensemble p = exp[(3{H — fJ-N)] where 
N is the total-particle-number operator, the fluctuation- 
dissipation theorem states that F and p are related by 

Faa{Xo;LL>,p) = -i[l/2 - nF-D{uj - p)]paa{Xo;i^,p). (3) 

In Fig. m we present, at the late time Xq = 18.9 n^^m, 
the fraction / = iFf^{Xo;ijJ,p)/ p^f{Xo;io,p) as a func- 
tion of the frequency ui, for five different momentum 
modes p^. The lower left panel shows the inverse-slope 
function ln[(l/2 — f)~^ — 1] of / for run C, the lower 
right panel for run D. / is shown in a region including 
the spectral peaks where the argument of the logarithm 
is positive. Outside this region, it oscillates around zero, 
owing to the finite total evolution time after the quench. 
In run G, we find, as in B (not shown), that over the 
region of relevant uj this function is a straight line and 
therefore reflects a Fermi-Dirac function. Hence, accord- 
ing to the fluctuation-dissipation theorem, the system is 
thermalized over the depicted range of energies, despite 
the power-law tail in run G found for Wp > 30 wf, sec 
Fig.lSK. This can be understood by considering the spec- 
tral function in Fig. 3] (upper row) on a logarithmic scale. 
A second peak at negative frequencies, not present in an 
ideal gas, picks up extra contributions from the Fermi 
sea. It is this second peak which causes the power-law 
overpopulation at high momenta, reminiscent of the BGS 
zero-temperature depletion of the Fermi sea in a weakly 
interacting system. Hence, from run G one may precon- 
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